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FOREWORD 


This report describes work performed by the Pratt & Whitney Aircraft 
Group, Commercial Products Division, United Technologies Corpora- 
tion, for the National Aeronautics and Space Administration, Langley 
Research Center under Contract NAS1-14462. During this effort Mr. 
Lawrence E. Putnam was the NASA Project Manager and Dr. W. M. 
Presz, Jr. was the P&WA Program Manager. The report includes a 
user’s manual for an afterbody viscous flow computer program devel- 
oped under the contract. The user’s manual is presented as an appen- 
dix. 
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AN IMPROVED ANALYTICAL MODEL OF THE SEPARATED 
REGION ON BOATTAIL NOZZLES AT SUBSONIC SPEEDS 


by 

Walter M. Presz Jr., Ronald W. King and John D. Buteau 
Pratt & Whitney Aircraft, East Hartford, CT 06108 
SUMMARY 

Boundary layer flow separation often occurs on the afterbody of aircraft, increasing drag and 
decreasing performance. Although the separation phenomenon is theoretically understood, 
the complexities of the governing equations limit solutions to simplified flowfields. Because 
present methods of predicting and analyzing afterbody flow separation are not reliable, the 
current program was undertaken. 

A practical engineering calculation has been developed to model the viscous effects of a sep- 
arated, reverse-flow region on afterbody pressures and drag. This viscous calculation is de- 
signed to be iteratively coupled with any existing inviscid-flow calculation through an aero- 
dynamic interface. A standard boundary layer displacement thickness is used to modify the 
afterbody shape where the flow is attached. A discriminating streamline calculation is de- 
veloped to account for the displacement effects of the reverse flow in separated regions. The 
discriminating streamline calculation includes the effects of skin friction, axial pressure grad- 
ients, and jet entrainment. Skin friction was found to have little effect on the shape of the 
discriminating streamline. Axial pressure gradients and jet entrainment effects, however, do 
drastically change the shape of the discriminating streamline from the conical shape pre- 
viously used. 

The viscous flow calculation was coupled to a potential flow calculation. Predictions ob- 
tained with this computer program were compared with experimental data for afterbodies 
with solid jet simulators previously tested by NASA. The analysis was found to accurately 
predict both the magnitude of the measured afterbody pressures and drag and also their 
variation with Reynolds number, Mach number, and afterbody shape. Predictions for cases 
with a flowing jet were also compared with a limited amount of available data. The code 
accurately predicts a double cell reverse flow region near the afterbody-jet juncture, and 
reasonably predicts the afterbody pressures. Recommendations for new data required to 
improve the accuracy of the technique predicting the effects of the jet on afterbody pres- 
sures and drag are also presented. 


INTRODUCTION 

Boundary layer separation frequently occurs on the afterbody of jet aircraft, increasing drag 
and decreasing engine performance. Currently employed methods of predicting and analyzing 
this phenomenon are not reliable, which handicaps the designer and increases the uncer- 
tainties in predicting performance. To ameliorate this condition, this research program was 
undertaken. The program is based on earlier work conducted by Dr. W. M. Presz of Pratt & 
Whitney Aircraft and Dr. E. T. Pitkin of the University of Connecticut (ref. 1 & 2). 



Boundary layer separation, a fundamental problem of aerodynamics, is caused primarily by 
two related factors: viscosity and an adverse pressure gradient. Viscosity depletes the mo- 
mentum of the boundary layer flow in an adverse pressure gradient. At some point, because 
of the adverse pressure gradient, the flow breaks away from the wall, and a local area of 
reversed flow is formed. The point along the wall dividing the forward and reverse flows is 
defined as the separation point. 

An exact flow solution for boundary layer separation can theoretically be obtained from the 
finite difference calculations of the fully turbulent, compressible Navier Stokes equations. 
The required calculations, however, have only been solved for simplified flowfield cases, 
and have required such long computer times as to render the approach unrealistic at this 
time for engineering requirements. 

Current engineering techniques for predicting the separation point generally employ only 
attached boundary layer and inviscid-flow pressure calculations. A usual assumption of 
these techniques is that both the momentum and pressure gradients normal to the solid 
boundary contour are negligible in a thin flow region adjacent to the contour. This approxi- 
mation allows the momentum equations in a direction normal to the contour and several 
related flow-direction terms of the Navier Stokes equations to be omitted, resulting in the 
boundary layer momentum equation. The momentum integral equation is simply the inte- 
gral of the momentum equation through the boundary layer in the normal direction, and is 
limited, therefore, to the same basic assumptions. When either of these equations is used to 
calculate boundary layer development on a curved wall, the development is actually for a 
flat plate with the characteristic pressure distribution of the curved surface superimposed. 

Neither the differential boundary layer equations (finite difference) nor the momentum inte- 
gral equations have been very successful in predicting flow separation (ref. 3). A local singu- 
larity in the calculation occurs near the separation point, rendering finite difference tech- 
niques in gross error. In addition, all available turbulence models are based on attached boun- 
dary layers and become very inaccurate near separation. 

The momentum integral equations are inadequate because they require the calculation of a 
small difference of large numbers as separation is approached. Furthermore, the empirical 
skin-friction model used in most turbulent momentum integral calculations never reaches 
zero — and on a smooth axisymmetric surface with a large radius of curvature, the separation 
point is always associated with a zero wall shearing stress. In addition, the velocity profiles 
used in the turbulent momentum integral calculations are based on attached boundary layer 
assumptions, which are unrealistic near separation. 

Calculations have shown, however, that with proper assumptions (viz, a priori displacement 
thickness distribution through separation) the singularity can be eliminated in the finite dif- 
ference approach and the boundary layer calculations will yield valid approximations to the 
full Navier Stokes equations. The boundary layer assumptions are not, therefore, the critical 
limitations for predicting separation. Current techniques for calculating separation have 
been unsuccessful because separation and the resultant recirculation region drastically alter 
the afterbody pressure distribution, which in turn significantly alter the position of the sepa- 
ration point. Standard boundary layer techniques that use only the attached boundary layer 
and inviscid flowfield calculations to predict separation cannot account for this effect. 


2 . 
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The engineering approach for predicting flow separation on afterbodies developed during 
the earlier work (ref. 2) is iterative and combines all four major flowfield components shown 
in figure 1 : the inviscid flowfield, the attached boundary layer, the point of separation, 
and the separated recirculating-flow region. The iterative scheme first assumes that there is 
no separation. Based on this assumption, the inviscid flowfield and corresponding pressure 
distribution on the afterbody are determined. Next the boundary layer development is cal- 
culated by a turbulent boundary layer program. The separation point is then determined by 
a control volume calculation; and a recirculating-flow region is introduced downstream of 
the separation point. 

The streamline separating the external and recirculating flows is called the discriminating 
streamline and is treated as if it were the boundary of a body extending between the separa- 
tion point and the reattachment point on the real body. The discriminating streamline in re- 
ference 2 is assumed to be straight, and this assumption is used to define the reattachment 
point. A recirculating-flow velocity profile is assumed and combined with a mass and a 
momentum balance for the recirculating region to determine the angle this streamline makes 
with the afterbody at the point of separation. 

Once the separation and reattachment points have been determined, the inviscid flowfield is 
recomputed for the flow over the quasi surface as defined by body contour plus the discrim- 
mating streamline. The entire process is now repeated, and the reiterations continued until 
the solution converges. 

Although standard techniques were shown to be capable of accurately calculating the attached 
boundary layer and the inviscid flowfield, new flow models had to be developed to calculate 
the point of separation and the recirculating flow. The earlier work (ref. 2) showed that 
blockage is the dominant effect of a separation region on an afterbody flowfield. Once the 
point of separation and a discriminating streamline designating the displacement effect of 
the recirculating region are determined, afterbody pressures and drags can be calculated from 
an inviscid analysis of an aerodynamic contour defined by the body and the viscous flow dis- 
placement thickness of both the attached boundary layer and the recirculation region. 

Available experimental pressure data typically indicated a constant pressure region covering 
the majority of the model-afterbody area between the separation and reattachment points. 

It was determined from both analysis and data that a conical discriminating streamline would 
produce the observed constant-pressure region. A conical, discriminating-streamline model 
was utilized in reference 2. 

Although agreement is good between the analyses and data presented in references 1 and 2, 
available data is limited. Furthermore, the flow calculation developed has several inherent 
problems in the flow which could limit its application: the discriminating streamline shape 
is crude and is based on insufficient data; the no-slip wall boundary condition is not accounted 
for in the recirculation region; and jet entrainment effects are neglected. As a result skin 
friction forces cannot be calculated through the recirculation region and drag could be inac- 
curate on an afterbody with a flowing jet. The program described in this report addresses 
these limitations, and identifies new data required to develop a better flow model for after- 
body flow separation. 
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SYMBOLS 
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P 

P t 

q 

R 

R 

e 

R 


r. 

j 


s 

s x 

T 


u 


u. 

J 


V 

( 

T 


V. 

J 


cross-sectional area 

minimum cross-sectional area of afterbody 
maximum cross-sectional area of afterbody 
equivalent length of fiat plate boundary layer 
constant representing the wake spread rate 
afterbody pressure drag coefficient: ^ / C dA 

skin friction coefficient m 

static pressure coefficient, (P - P^j/q 

maximum diameter of model 
base diameter of afterbody 
sting diameter 

boundary layer shape factor, 5*/# m 

height of discriminating streamline above body surface in the s coordinate system 
axial momentum in the X direction for axisymmetric flow 
axial momentum in the X direction for two-dimensional flow 
length of model from nose to start of afterbody 

length of afterbody 
width of wake 
width of shear layer 
Mach number 

normal distance from afterbody surface 
local static pressure 
local total pressures 
freestream dynamic pressure 

Reynolds number based on distance from noise to start of afterbody 

local Reynolds number based on axial distance from model nose 

local Reynolds number based on boundary layer momentum thickness 

radial coordinate system 

radius of outer shear layer edge 

radius of the discriminating streamline 

radius of outer edge of wake 

radius of inner shear layer edge 

coordinate system perpendicular to streamline in the recirculating region 

surface area defined by the separated-flow region 

static temperature 

total temperature 

local axial velocity 

maximum reverse flow velocity 

inviscid jet velocity 

local inviscid velocity 

radial velocity at edge of shear 

radial velocity at edge of w ake 
radial velocity of jet 
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SYMBOLS (Cont'd) 

X 0 - coordinate system parallel to flow streamlines in recirculating region 
x - axial coordinating system 

w - local recirculating velocity 

a - angle at profile inflection point 

a - angle between the afterbody surface and the profile inflection point 

W 

j3 - afterbody mean angle 

v - kinetic viscosity 

5 - boundary layer thickness 

5* - boundary layer displacement thickness 

7 - nondimensional coordinate, (7 c -7)/® 2 

\jj - angle of the discriminating streamline 

p - density 

o - jet spreading rate parameter 

r - local shear stress 

6 - local flow angle 

0 - boundary layer momentum thickness 

m 

f nondimentional coordinate, (7 - -y c )/C 1 

"n - unit normal vector 

s - vector having magnitude equal to the control volume surface area and a direction 
parallel to the body surface 

Subscripts: 

1,2 - station indicators 

av - average parameter 

body- along contour 

c - outer edge of jet shear layer 

d - along discriminating streamline 

e - local inviscid flow properties 

j - inner edge of jet shear layer 

n - point of minimum pressure on an afterbody 

o - start of the afterbody 

s - condition at afterbody separation point 

00 . freestream condition 


PROGRAM DESCRIPTION 

This research program conducted under NASA contract NAS 1 - 1 4462 , entitled “An After- 
body Viscous Flow Program”, has extended the usefulness and increased the accuracy of 
the analytical procedure developed earlier by Dr. Presz and Dr. Pitkin. The analysis itera- 
tively combines four major components: the inviscid flowfield, the attached boundary lay- 
er, the point of separation, and the separated recirculating flow region. The work described 
herein improves the recirculating flow model to more accurately predict the effects of sepa- 
ration on axisymmetric nozzle afterbody performance. Previously this analysis had not ac- 
counted for the effects of varying density, axial pressure gradients, skin friction, and jet en- 
trainment. These effects are now included. 
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The technical effort was organized into four tasks. During Task I the straight-line, discrimin- 
ating-streamline analysis of reference 2 was further developed to include the effects of vary- 
ing density, pressure gradients, and skin friction. This was accomplished by reformulating 
the analysis procedure to a discrete control-volume calculation. During Task II an entrain- 
ment model was added to the procedure, accounting for jet effects on separation. A viscous- 
flow computer program was developed during Task III by combining available, attached- 
boundary-layer and separation-point calculations with the new, discriminating-streamline 
model. The accuracy of the procedure was assessed during Task IV. 

TASK I - EXTENSION OF METHOD TO ACCOUNT FOR SKIN FRICTION EFFECTS 

Although the conical discriminating model developed by Dr. Presz and Dr. Pitkin provides a 
method for realistically predicting separated afterbody pressures and drag (ref. 1 & 2), the 
procedure has several limitations: separation angles predicted at Mach numbers above 0.8 
are too large (see figure 2), and the predicted recompression pressures are often much too 
high near the afterbody-sting junction (see figure 3). Both of these errors are believed to be 
due to the effects on the discriminating-streamline shape of variable density, axial pressure 
gradients, and wall skin friction. The objective of Task I was to correct the deficiencies by 
developing a discrete-control-volume analysis for the wake which would include these effects. 

A conceptual model of the flowfield in the vicinity of a separated afterbody is shown in 
figure 4 along with the discrete control volume used to calculate the discriminating-stream- 
line shape. The upper surface of the control volume is inclined at an angle \p with the wall 
and, since the upper surface is a flow streamline, no mass flux passes through it. The recir- 
culating-flow streamlines are assumed to be source-like in direction, with each streamline 
having a common origin. This assumption allows a smooth transition from the flow direc- 
tion on the body to that on the discriminating streamline. The inside walls of the calculation 
control volume are constructed normal to the streamlines. 


Applying the continuity and momentum equations to the axisymmetric control volume, one 
obtains 


f 1 

/ (p w )i. 

J n 


27 r / (pw) t (rj +s) ds - 27 r J (p w) 2 ( r 2 + s ) ds = 0 

o 


/' 

j o 


eq. (1) 


and 


rd? x ~P B ( 27r r av, ) AX o + '(r d ) 2n ( 


where 


•/' 


* 


J x = 27 t j p w 2 (cos — s ) (r + s) ds 


f av + h av ) COS ^av AX 0 = J Xi2 -J X)1 e <l-( 2 ) 


eq. (3) 
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and 


h ^ 2 

r 


r = 


^ 1.0 


r, + r~ 


2 


- tl 111 

= 2 

^X o 0 

= X 0 4 , 

= (hj +h 2 )/2 

The recirculating flow region can be considered a free-shear layer. Nash (ref. 4); and Vasilu 
(ref. 5) have used the similarity between the velocity distribution of a free-shear layer and 
the velocity of a jet flow in the analysis of problems involving shear layers associated with 
boundary layer separation. The free-shear layers encountered here are treated in like man- 
ner. Thus, the velocity across the reverse flow regions is assumed to vary as an error func- 
tion. The exact profile used is a modification of one developed by Korst and Chow (ref. 6) 
and which has been shown to accurately model isoenergetic turbulent mixing between two 
compressible streams at constant pressure. 


* 


av 


S 

h 

h 


eq. (4) 

where 


a = 12 (1 +.23 M e ) 



Equation (4) represents the recirculating-flow velocity as a function of streamline coordin- 
ates s and X 0 . This equation provides a continuous variation in both magnitude and direc- 
tion for the velocity between the reverse flow at the wall and the inviscid stream velocity. The 
angle 9 represents the velocity direction measured with reference to the local body surface. 

The velocity profile is shown in figure 5 ; the maximum velocity in the reverse-flow direction 
is u B . A laminar body layer is superimposed on the reverse-flow region, and is assumed to 
grow from the trailing edge of the afterbody for a flowing jet and from the reattachment 
point for a solid sting. The exact solution for laminar boundary layer growth on a flat plate 
is used for the boundary layer profile, as indicated in figure 5. Again, the angle inclination 
between the discriminating streamline and afterbody surface is i //. It is clear from the geom- 
etry in the figure that 

\p = aj +a w eq. (5) 


w = ^ + ^ ^ 1 + erf [a (0 - <* w ) ]'^ u B 
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The following relationships for the turbulent shearing stress along the discriminating stream- 
line and the angle a, can be derived from the tables presented by Korst and Chow (ref. 6). 

a 

eq. (6) 


U T> U„ 

0.1481 -0.0478 M -0.1278-^+ 0.1632 M +0.3990 

e U e e U p 


'd 

P~A " 




and 

0.2090 + 0.0226 M, +0.3080 

a = e eq. (7) 

a 

Equations 1,2, and 3 and the remaining defining equations comprise a proper set of three 
equations with three unknowns: h, a w , and u fi . Given the point of separation, body shape, 
and inviscid flow thermodynamic properties in the axial direction, the shape of the discrimin- 
ating streamline can be obtained. The equation set includes the effects of compressibility, 
axial pressure gradient, and skin friction. 

Figures 6 through 14 were obtained by applying this analysis to nine test cases chosen 
jointly by NASA and P&WA. These figures show the calculated pressure and separation loca- 
tion of separated flow regions for three afterbody contours at four Mach numbers. This data 
is from pressure and separation location experiments by Reubush (ref. 7) and Abeyounis 
(ref. 8), respectively. The results shown in the figures were obtained with an axisymmetric 
cone-cylinder forebody (L/D = 8) tested with various circular-arc afterbodies with solid- 
plume simulators, the Three afterbodies analyzed have mean angles and 2/D ratios of 17.0° 
and 0.8, 13.8° and 1.0, and 1 1.0° and 1.0. For these calculations the separation location 
was set equal to the experimental values obtained by Abeyounis and the pressure distribu- 
tions from experiments by Reubush. In these figures the straight-line predictions based on 
the reference 2 analysis are also shown for comparison. The discrete control volume results 
show considerable curvature in the shape of the discriminating streamline and significant 
variation with separation Mach number. Equivalent separation angles were obtained for the 
discrete-control-volume results by drawing a conical surface between the point of separation 
and the predicted flow reattachment point. A comparison of these equivalent separation 
angles and the reference 9 data is shown in figure 15. As can be seen, the results obtained 
with the discrete control volume analysis are in better agreement with the data than are the 
results obtained with the reference 2 analysis. 

Figure 16 shows the effects of the skin-friction, axisymmetric-flow, and pressure terms from 
equations 1 through 7 on the shape of the discriminating-streamline. Four distinguishably 
different discriminating streamlines are shown: (1) the straight conical streamline determined 
by the method of reference 2; (2) the predicted streamline, using equations 1 through 7 but 
neglecting skin-friction, pressure, and axisymmetric-flow terms; (3) the predicted streamline 
using equations 1 through 7 but neglecting only skin friction and pressure terms; and (4) the 
predicted discriminating streamline, using equations 1 through 7 both with and without skin- 
friction terms. 
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The predicted discriminating streamlines in figure 1 6 are essentially the same with a C f - 
0.003 and with a C f = 0.000. This is consistent with equation 2 where skin-friction and 
shear-stress term order-of-magnitudes can be compared. 

We have from equation 2 

C 

skin friction term Pb u b 2 ^ 7r r av) 

shear stress term r d 2 tt (r ay + h av ) cos \p av 


Since 


p u ‘ 
e 


and 


Therefore 


r av + h av ” r av 


COS i// av « 1.0 

Pe * Pb 


skin friction term 
shear stress term 


(x) 


< < 1.0 


Thus, skin friction effects on the discriminating-streamline shape should be second order, 
which agrees with results shown on figure 16. Figure 16 also shows, however, that the effect 
of axisymmetric flow tends to increase the separation angle and, thus, the size of the separa- 
tion bubble. The pressure-gradient term increases the discriminating-streamline angle in a 
favorable gradient region and decreases the angle in an adverse gradient region. 

A potential flow analysis (ref. 10) was used to calculate the flow over the various discriminat- 
ing streamlines presented in figure 16. Figure 17 compares the predicted pressures for con- 
figuration 1 as functions of axial location. These results show that the discrete-control-vol- 
ume analysis can accurately predict the effects of separation on afterbody pressures. 

Figure 1 6 shows that even when the skin-friction, axisymmetric, and pressure terms are 
dropped there is still a considerable difference between the discriminating streamline obtained 
from equations 1 through 7 and the one presented in reference 2, which was obtained using 
similar equations but with constant density and a constant separation angle assumed through- 
out the recirculation flow region. This discrepancy is investigated in the following paragraphs. 
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By assuming constant density, letting r get very large compared with h (i.e., to approach two- 
dimensional flow), and dropping the pressure-gradient and skin-friction terms, equations 1 
through 7 reduce to 


,,, f 1 (pu), ao - x„,2 [ 

J n J n 




(pu), d0 = 0 


eq. (8) 


0.1481 - 0.0478 (M ) - 


0 1278 (x) + 0 1632 Me (^) av + 0 ' 3990 (^f) a 
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eq. (9) 


eq. (10) 


where 


K. = X 



p u' 2 (cos d) dd 


X 


O 


h_ 

4> 


eq. (11) 


Equations 8 through 1 1 comprise a proper set of equations for calculating the angle of a 
conical surface rather than a discrete discriminating separation streamline, and are similar 
to those presented in reference 2. Figure 18 provides a comparison between the results 
obtained with the simplified equation set which assumes incompressible flow and small- 
angles, as in reference 2, and those obtained by solving exactly the equations using numeri- 
cal techniques. The results show that both compressibility effects and errors associated with 
small-angle assumptions can have a large effect on the predicted discriminating streamline 
angle. 
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The complete, exact solutions of equations 1 through 7 have been programmed and coded 
for one part of the complete separation analysis shown in figure 19. For the remaining com- 
ponents, existing flow calculations are used. A user’s manual for the computer program de- 
veloped is presented in the appendix. 

TASK II - INCLUSION OF JET ENTRAINMENT EFFECTS 

During Task II, the analytical model developed during Task I was extended to include jet en- 
trainment effects, which are neglected by currently utilized models. The normal procedure 
has been to remodel the jet by using a solid sting, a procedure that can lead to large drag 
errors (ref. 1 1 ). 

Two approaches for analyzing the separation bubble-jet shear-layer region were examined dur- 
ing Task II. The first and simpler of the approaches was to determine the magnitude of the 
entrainment effects on the calculation procedure by imposing a slip-velocity and a shearing 
term at the trailing edge of the afterbody. (See figure 20). This method can be expected to 
express the jet effects quite accurately near the afterbody trailing edge. However, as the jet 
mixing layer grows and affects the profile of the recirculating region, the method becomes 
increasingly in error. As shown in figure 21 , the effect of the jet shearing term on afterbody 
pressure is quite large. Thus, neglecting the jet mixing layer cannot be justified. For this 
reason the second approach was developed. The second approach is more complete and in- 
cludes the jet mixing layer. 

The model for the second approach assumes an incompressible, turbulent boundary layer and 
that the radial pressure gradients are negligible. By satisfying integral forms of the equations 
of motion, the developed procedure calculates shape parameters for an assumed mean-velocity- 
profile. 

The discriminating streamline is then calculated and used as input for the inviscid-flowfield 
calculation procedure. Communication of jet entrainment to the separated flow over the 
nozzle afterbody is assumed to take place by means of changes in the inviscid flowfield that 
result because of the effects of jet entrainment on the discriminating streamline. 

The schematic shown in figure 22 represents the flowfield assumed for a nozzle afterbody 
with both separated flow and a high-pressure jet. The separation bubble exists beyond the 
nozzle exit plane and consists of two recirculation zones. One recirculation zone is associa- 
ted with the freestream flow; the other with the jet-stream flow. 


A schematic of the flow model used to simulate the flowfield of the jet separation bubble is 
shown in figure 23. The viscous region is assumed to consist of two parts: a jet shear layer 
and a separation bubble wake. Velocities at the outer and inner edges of the viscous region 
are assumed to be equal to the inviscid freestream velocity and to the inviscid jet velocity, re- 
spectively. The minimum velocity of the profile is assumed to be at the interface of the jet- 
shear layer and the separation wake. The discriminating streamline for the freestream flow 
must exist between the outer edge and inner radius of the wake region. 
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The velocity profile selected to simulate the viscous region is shown in figure 24. The profile 
consists of two turbulent mixing layers which treat the minimum velocity as a counter flowing 
stream. The expressions describing the velocity profiles are modifications of the polynomial 
function used extensively by Abramovich (ref. 12). 

For the jet shear layer — Oj < r < r c ) 


u = - u B + A Uj (2 tj 1,5 — t? 3 ) 


eq. (12) 


U j + U B 


For the wake region — (r,<r<r e ) 

U = -u B + Au e (2f us -f 3 ) eq. ( 1 3) 

where Au e = u g + u B 

f = r — r n 

J C 


One additional geometric parameter is shown in figure 24: the radius of the discriminating 
streamline (r d ). 

Two of the governing equations for the jet entrainment model were obtained by integrating 
the differential forms of the equations of motion and applying the Leibniz formula to the 
integrals. 

The resulting continuity equation for the viscous flowfield is 


where 


Au- = 


9 

9x 




e dr j dr e 

Urdr+U i r j^- U er eir 


+ V e r e “ V j r j = 0 


eq. (14) 


Assuming that the incompressible form of Bernoulli’s equation is applicable results in the 
following expression for the axial momentum equation 


9 

9x 


/ r e d^ 2 

u2rC * r +U j 2r j - U J V j r j - U e r . 

1. 


dr 


e 

e dx 


r e du e 

+ u e v e r e u e — =0 eq. (15) 
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Two additional equations were added to the equations of motion to complete the set of 
governing equations. One equation is the continuity equation for the shear layer 



where v c 

The other equation 


drj 

urdr + u.r. — 
J J dx 


- ur 
m c 


dr 

c 

dx 


+ Vc'Yj = 0 


— radial velocity at the outer edge of the shear layer 
is the continuity equation for the discriminating streamline 


eq. (16) 


d 

dx 



ur dr + UjTj 



0 


eq. (17) 


Using the assumed velocity profile with the governing equations results in the following set 
of eight unknowns 


dUg 

dx 

d ± 

dx 

d£ 2 

dx 

dij 

dx 

dx 


axial rate of change of the minimum velocity 
slope of inner edge of shear layer 

spreading rate of the shear layer 

spreading rate of the wake 

slope of the discriminating streamline 

radial velocities at points within the viscous region. 


Since there are eight unknowns and only four equations, four additional conditions have to 
be specified to obtain a solvable set of equations. 

Since the jet shear layer is assumed to dominate the viscous flowfield, three of the necessary 
closure conditions relate to the characteristics of the shear layer. The first assumption is that 
the shear layer spreads at a constant rate, as suggested by Abramovich (ref. 12). 

d£, 

=0.25 eq. (18) 

dx 
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The second assumption is that the radial velocity at the inner edge of the shear layer is zero. 


Vj = 0 eq. (19) 

The third closure condition is that the net flow through the outer edge of the shear layer is 
defined by 


dr c 

- u b r c + Vc = k: Vc 


eq. (20) 


where k is a constant whose value is selected to maintain a continuous discriminating stream- 
line slope at the initial jet entrainment station. 

The final closure condition assumes that the outer wake region has a constant spreading rate. 
The value of the spreading rate is assumed to be given by 


where 



dx 

dx 


dx 


+ 0.25 = C 
e 


= slope of the nozzle afterbody at the nozzle exit, 
e 


eq. (21) 


TASK III - DEVELOPMENT OF A VISCOUS FLOW COMPUTER PROGRAM 

A computer subroutine package was developed during Task III. This package, called VISCUS, 
calculates boundary layer flow separation by utilizing the procedures previously described 
and by interfacing them with existing programs that calculate the external inviscid flow. The 
subroutine package is described in the appendix. 

To define the equivalent aerodynamic contours of afterbodies with separated flow, the dis- 
criminating-streamline wake calculations of Tasks I and II were combined with: 

1 ) a separation point prediction method 

2) a modified Reshotko-Tucker boundary layer analysis (ref. 13). 

3) a displacement thickness, inviscid flow, interface calculation for attached boun- 
dary layers. 

The calculation procedure, shown schematically in figure 25, was developed for iterative 
matching with an existing inviscid flow procedure. The viscous flow procedure calculates 
the aerodynamic contour based on a given flow-pressure distribution, and the inviscid pro- 
cedure calculates the pressure distribution based on the aerodynamic contour. The iteration 
between the two calculations continues until the aerodynamic contour and the pressure dis- 
tribution are in equilibrium. 
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The basic assumptions in the analysis are: 1) the inviscid or external flowfield imposes its 
axial pressure field on the boundary layer and separation region, and 2) the normal pressure 
gradients in the boundary layer and separation region are not important. The following 
paragraphs describe each of the major components of the VISCUS subroutine. 

Separation Point Prediction 

The VISCUS subroutine package provides the “user” with four options for calculating the 
separation: Stratford’s criterion. Page’s criterion, Goldschmied’s criterion, and the control 
volume criterion. 

Stratford’s Criterion (ref. 14) — This approach is based on an approximate solution of the 
equations of motion; a single empirical factor is required. The equations. are integrated using 
a modified inner and outer solution technique. It is assumed that the outer part of the boun- 
dary layer is affected only by the initial velocity profile and the downstream pressure gra- 
dient, and the inner part of the boundary layer is locally in equilibrium and is independent 
of upstream conditions. At a Reynolds number of the order of 10 6 , the criterion is: 


C p ,s fbd(C p >s )/db] % = 0 .98(R eb )°l e Q- (21) 

The distance b represents an equivalent length of flatplate, constant-pressure boundary-layer 
growth. Re^ is the Reynolds number based on this equivalent length. 

Page’s Criterion (ref. 15) — This is probably the most commonly used approach for pre- 
diction of the local separation point. The recirculation region is modeled as a base-type region 
and entrainment equations for constant pressure viscous mixing are used to derive an equili- 
brium pressure. In its original form for subsonic flow, the criterion C s = 0.38 where C s is a 
modified separation pressure coefficient defined as 

C s = (2(P S -P 00 )/ 7 P S M S 2 ) eq. (22) 


Eilers (ref. 1 6) modified this criterion for use on axisymmetric afterbodies where the boun- 
dary layer goes through a strong expansion before recompressing to a separated region. This 
modified version, which will be used herein, is C g = C n + 0.38, where C n is defined as 

C n = (2(P n - PqqI/t P n M n 2 ) eq. (23) 


The subscript n indicates reference to the point of minimum pressure on the afterbody. 
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Goldschmied’s Criterion (ref. 1 7) — This criterion was originally derived for planar incom- 
pressible flow. The criterion is based on a line of constant total pressure occurring in a boun- 
dary layer, even under an adverse pressure gradient. At separation the height of the constant 
pressure line and the laminar sublayer are equal. The resulting criterion becomes C p> s = 

200 C fh where C f n is the skin friction coefficient at a minimum pressure point on the 
afterbody. This separation prediction technique is independent of the pressure distribution 
downstream of the minimum pressure point and is dependent on Reynolds number and Mach 
number only through the skin friction coefficient. 

Control Volume Criterion (ref. 18) — A control volume in the boundary layer between the 
minimum pressure point and the point of separation on the afterbody is used. Conservation 
of mass and momentum through the volume, along with the compressible law of the wall 
law of the wake boundary-layer profiles are used to calculate the separation pressure. 


Boundary Layer Analysis 


The attached, turbulent boundary layer is calculated using an integral approach. The need 
to assume something explicitly about the local Reynolds stresses is eliminated by integrating 
the flow equations across the boundary layer. In this manner, the following momentum inte- 
gral equation is derived. 


dd_ 


dx 




(2 + H) + 




eq. (23) 


A solution of this momentum integral equation for turbulent boundary layers is obtained 
using the approach and equations of Reshotko and Tucker (ref. 13). An approximate shear 
integral obtained by a correlation of extensive data in an adverse pressure gradient is em- 
ployed (ref. 19). 


Aerodynamic Interface 


An aerodynamic contour is calculated by adding a displacement thickness to the basic con- 
tour where the boundary layer is attached and to the discriminating streamline where the 
flow is separated on the afterbody. The contour is then smoothed between each iteration. 


The smoothing routine uses finite difference equations with a Taylor series expansion to de- 
termine, among other properties, the curvature at the midpoint of each five point interval 
describing the body. .Curvature is calculated both with and without the midpoint. When the 
curvature is larger with the midpoint not included, the curve is under-constrained with the 
four points. Smoothing or moving the centerpoint would significantly alter the local char- 
acteristics of the curve near the midpoint and, therefore, no smoothing occurs. If the curva- 
ture is smaller with the midpoint, the midpoint is adjusted in a systematic manner to increase 
the curvature. Thus, a numerical model with a smooth, continuous, second derivative is as- 
sured. All sharp comers or slope discontinuities are faired out. 
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TASK IV - ASSESSMENT OF ACCURACY OF METHOD 

The afterbody viscous flow calculation from Task III was interfaced with an existing, inviscid, 
potential-flow program that uses a combination of the Neumann method, (ref. 10) and the Lab- 
rujere compressibility correction (ref. 20).i The complete analysis was applied to several after- 
body contours and flow conditions, and the results compared with available data to deter- 
mine the accuracy of the calculation methods. The test cases included fifteen chosen jointly 
by NASA and Pratt & Whitney Aircraft. The comparisons, discussed below, have been di- 
vided into afterbody flowfields with solid-plume simulators and afterbody flows with a jet: 
the solid-t»lume comparisons relate to Task I calculations; the flows with jets, to Task II cal- 
culations. 

Solid Plume Results 

A jet exhaust plume affects afterbody boattail pressures in two ways: through blockage and 
through entrainment. A significant amount of data for separated afterbody contours has 
been recently obtained by using solid jet simulators, or rearward stings (ref. 7,8,9, 11,21 
and 22). Although this type of jet simulator accounts for only the blockage effect of a jet, it 
does provide a useful media for studying the separation region over an afterbody, including 
the reattachment zone on the sting. A study of the accuracy of the analytical prediction 
methods developed in Task I was determined using this data. This study included variations 
in afterbody pressures and drag as a function of Reynolds number, Mach number, afterbody 
contour, and separation location. Results obtained by varying each of these parameters will 
be discussed in the following paragraphs. 

A thorough experimental study of afterbody pressures and drag as a function of Reynolds 
number has been conducted by NASA (ref. 21) in the NASA-Langley 1 /3-meter cryogenic 
tunnel. A comparison of the analytically predicted pressures with data from these tests are 
shown in figure 26 for a circular arc-conic afterbody at a Mach number of 0.601 and a 
Reynolds number of 4.29 x 10 7 . The upper part of the figure presents the afterbody con- 
tour. Several iterations of the pressures predicted by the analysis are presented. In the first 
iteration, the calculation procedure solves for pressures over the actual contour without a 
boundary layer. A boundary layer displacement thickness is used to modify the afterbody 
contour on all subsequent iterations. In the third iteration, the procedure solves for the 
pressures over the afterbody contour modified to also include a constant-pressure discrim- 
inating separation streamline. In the fifth and subsequent iterations, the procedure solves 
for the pressure over the afterbody with the discrete discriminating streamline calculated to 
include the effects of skin friction and pressure gradients. The separation point is averaged 
with the previous value, and the thickness of the boundary-layer displacement is averaged 
with the previous two values after the fifth iteration. (See appendix for details of this itera- 
tion procedure.) The convergence shown in figure 26 is excellent after seven iterations for 
both predicted pressure distribution and separation-point location. The separation point 
results produced in figure 26 are based on the control volume technique. Since the same re- 
sults were observed after the seventh iteration in other applications of the analysis, the iter- 
ation scheme was set to fix the separation point at the value predicted during the seventh 
iteration. This scheme accelerates convergence of the boundary layer, minimizes instabilities 
in calculations for more severely separated contours, and provides the best comparisons with 
available data. 
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A comparison of the results from the final (i.e., tenth) iteration with data for the same body 
is shown in figure 27. The agreement between analysis and data is excellent. The largest 
differences between predicted and measured pressures are over the initial portion-of the 
afterbody (x/D = 0 to 0.2); the reason for the discrepancy is not clear. The pressures in 
the highly viscous flow region at the trailing edge of the nozzle are matched quite well. A 
comparison of measured and predicted pressure drag (also shown in figure 27) shows that 
the prediction is higher than the measured pressure drag. This result has to be function of 
the pressure discrepancies over the initial portion of the afterbody. 

Figure 28 presents a comparison of measured and predicted pressures and drags over this 
same afterbody at the lower Reynolds number of 2.61 x 10 6 . The wide variation in Reynolds 
number in the data was obtained using the NASA Langley 1 /3-meter cryogenic tunnel. 

Again, agreement between analysis and data for the pressures is excellent except near the 
onset of the afterbody. The analysis predicts the weaker expansion and recompression on 
the afterbody, which occurs at this lower Reynolds number and also that the separation point 
moves downstream with increasing Reynolds number. The pressure drag is again over-estimated. 
In comparing figure 27 and 28, the analysis predicts a slight decrease in afterbody pressure 
drag with Reynolds number, and the data verifies this trend. 

Separated afterbody pressures have been measured on a full circular-arc afterbody (£/D = 

0.80 and d b /D = 0.51) at several freestream Mach numbers (ref. 7 and 11) and a comparison 
of the predictions with the data is shown in figures 29, 30, and 3 1 for Mach numbers of 
0.4, 0.6, and 0.8, respectively. The agreement between analysis and data is excellent at a 
Mach number of 0.4 and is good at 0.6 and 0.8. The analysis predicts that with an increase 
in Mach number there will be an increase in expansion, an upstream separation point move- 
ment, and an increase in drag; these predictions are all verified by the data. At the higher 
Mach numbers the recompression pressures are slightly overestimated by the analysis. The 
major reason for this discrepancy is in the separation-point calculation. This is emphasized 
in figure 32 where the pressures over the afterbody are again predicted by the analysis at a 
Mach number of 0.6, but with the separation point fixed at the experimental value. The 
agreement between analysis and data becomes much improved, substantiating that the sep- 
aration point calculation is the major sources of error for existing discrepancies between 
analysis and data. 

Additional comparisons were made for several different contours, and these results are pre- 
sented in figures 33 through 37. Figure 33 shows predicted results for a circular-arc after- 
body having an £/D =1.0 and a d^/D = 0.61 (ref. 24) at a Mach number of 0.8. Although 
the predicted separation location is slightly too far upstream, agreement between predicted 
and measured pressures is excellent. 

Figure 34 presents the results for a circular-arc, conical afterbody that was tested by 
Shrewsbury (ref. 23) and which has a ratio of radius-of-curvature to maximum-diameter of 
0.5 and a ratio of afterbody-length to maximum-diameter of 0.5715. This afterbody con- 
tour is very nearly conical, as can be seen in the upper portion of figure 34. The measured 
pressures near the juncture of the afterbody and the maximum diameter indicates a very 
rapid, localized expansion region. The analysis drastically under predicts this expansion 
region. Since this area is far upstream of any separation point or viscous region, this discrepancy 
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appears to be due to limitations of the inviscid flowfield calculation or the contour smoothing 
employed. Further downstream near the junction of the afterbody and sting, the analysis 
accurately predicts the pressures. 

A comparison of predicted and measured pressures is presented in figure 35 for a full cir- 
cular-arc afterbody having a L/D = 8,0, and £/D = 1.0, and a d b /D = 0.51 at a Mach number 
of 0.798— the contour is shown in the top of the figure. The agreement between both pre- 
dicted and measured pressures is again excellent. This figure also shows the effect on the 
predicted pressures as the separation point is artifically moved. When the separation point 
is moved downstream Ax/D = 0.05 from the experimental value, too much recompression 
near the sting-afterbody juncture and relatively little difference in afterbody drag results. 

When the separation point is moved upstream Ax/D = 0.05 from the experimental value, less 
recompression near the afterbody juncture and an increase in afterbody drag results. These 
results emphasize the importance of accurately predicting separation point location when us- 
ing a patched inviscid, viscous interaction procedure for estimating afterbody pressures and 
drag. 

Jet Plume Results 

During this portion of Task IV, the accuracy of the jet-entrainment model developed under 
Task II was assessed. Because detailed information is not currently available for jet/separa- 
tion bubble interaction flowfields, the accuracy had to be checked indirectly. This was 
accomplished by using the wake entrainment model to predict the effect of jet entrainment 
on afterbody pressures and comparing the results with data. The model was incorporated 
into the afterbody calculation analysis described previously. The procedure was then ap- 
plied to five of the test cases chosen jointly by NASA and Pratt & Whitney Aircraft. Ex- 
perimental data for the test cases— which included variations of nozzle pressure and free- 
stream Mach number— were obtained from an earlier NASA experimental study (ref. 1 1 and 
22 ). 

The iteration procedure used for analysis of the jet entrainment effects was identical to the 
one used for the solid-jet simulator analysis, except the jet entrainment model was used to 
calculate the discriminating streamline downstream of the nozzle exit. The inviscid core of 
the jet plume was assumed to be independent of the external flowfield, and the shape was 
assumed to be identical to that of an inviscid jet issuing into a quiescent atmosphere. 

The effect of jet entrainment on the discriminating streamline of a separation region is shown 
in figure 36. A comparison is presented for discriminating streamlines calculated for flow over 
a circular-arc afterbody at a Mach number of 0.6 and a nozzle total pressure ratio of 2.91. The 
dashed line is the calculated discriminating streamline, assuming the inviscid core to be a solid 
sting, thus accounting for plume blockage only. The solid line is the discriminating stream- 
line calculated with jet entrainment included. The separation point was not allowed to change 
for these calculations and was set equal to the experimentally observed separation point. A 
significant difference in the two streamlines can be seen at the end of the recirculating zone. 

In figure 37 the afterbody pressures calculated by the two modes are compared. These results 
indicate that the differences between the two discriminating streamlines downstream of the 
nozzle exit has a small effect on the afterbody pressures, implying that the mixing layer assumed 
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to grow along the discriminating streamline dampened out the effects of streamline variations 
near the end of the separated region. The results of the application of the analysis to the five 
test cases is presented in figures 38 through 42. The separation point was computed for these 
cases and was fixed at the value predicted in the seventh iteration. The shape of the discrim- 
inating streamline was also set as that calculated during the seventh iteration with subsequent 
iterations used to settle out the boundary layer superimposed upon the discriminating stream- 
line. The afterbody configuration analyzed for all of the test cases was a circular-arc after- 
body with 8/D = 0.80 and d^/D = 0.51. Figures 38 and 39 present comparisons of analyti- 
cal results with experimental data for a nozzle pressure ratio of 2.9 and Mach numbers of 
0.6 and 0.8, respectively. Figures 40 and 41 present a similar Mach number variation for a 
nominal nozzle pressure ratio of 2. Figure 42 presents the comparison for a Mach number of 
0.8 and a pressure ratio of 5. A good comparison between predicted and experimental separa- 
tion point can be observed for all cases. In all the comparisons the predicted measured pres- 
sure is good except right near the nozzle trailing edge. 

When a nozzle jet interacts with an external flowfield, the effective inviscid jet contour is 
changed (ref. 24) from what it would have been in a constant-pressure environment. Since 
all the results presented herein were based on an inviscid jet contour expanding into a constant- 
pressure environment, the calculations neglected suppression on the inviscid jet contour. And 
this could account for the small jet effects on the predicted afterbody pressures. Develop- 
ment of an iteration scheme that includes coupling with the calculation of the inviscid plume 
should be undertaken before further assessing the accuracy of the model. 

RECOMMENDATION FOR NEW DATA REQUIRED 

The primary area of uncertainty in the afterbody viscous flow calculation is the jet entrain- 
ment model. The analysis includes all the physics of the flowfield believed to be important, 
and yields a recirculation region that is consistant with available flow visualization and with 
intuition. Even so, pressures over the separated afterbody with a flowing jet represents the 
largest area of disagreement between calculation results and data. To resolve this problem, 
new data is required for improved calibration of the calculating procedure and for checking 
current assumptions made in the jet entrainment model. 

Calibration of the jet entrainment model requires test data that will enable refinement of 
the empirical-closure assumptions used to complete the set of governing equations. The 
specific closure assumptions that have to be refined are the spread and entrainment rates of 
the shear layer and the spread rate of the outer wake. To adequately define effective models 
for these flow characteristics, an elaborate test program is required. 

The program should be designed to define the separation-bubble, jet-entrainment region for 
axisymmetric nozzle afterbodies with high-pressure air jets. Variations of the flowfield with 
Reynolds number, Mach number, nozzle pressure ratio, and afterbody shape should be defined. 

The jet separation bubble interaction region should be mapped using optical techniques to pro- 
vide measurements of the pertinent flow characteristics. Use of conventional rake and probe 
or hot wire instrumentation is not recommended because the presence of the measurement 
device can significantly affect the flowfield in the recirculation region. 
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SUMMARY OF RESULTS AND RECOMMENDATIONS 

The conical discrimination streamline calculation for separated flows developed earlier by 
Presz and Pitkin (ref. 1) was extended to include effects of large separation angles, variable 
density, skin friction, and axial pressure gradient. Skin friction within the recirculating flow 
region was found to have a second order effect on the separated streamlines while pressure 
gradients were found to add a little curvature to the disciminating streamline shppe without 
drastically affecting its mean angle. The constant-density, small-angle assumption originally 
used in developing the conical discriminating streamline calculation was found to be a pri- 
mary reason for discrepancies between analysis and data at high subsonic Mach numbers. 

With these effects accounted for in the separation model, the analysis predicts the extent of 
the separation region on afterbody-sting models quite accurately for the data available. 

The discriminating streamline calculation was extended to also include jet entrainment ef- 
fects on the separation region. A wake-jet profile capability downstream of an afterbody was 
incorporated into the analysis. In this manner, jet entrainment was found to change the 
separation region on an afterbody from a single recirculating flow to a double counter-rotat- 
ing, recirculating flow. This results in the jet sucking down the discriminating streamline and 
drastically changing its shape downstream of the afterbody jet exit. 

The extended discriminating streamline calculation was combined with an attached-boundary- 
layer calculation, a separation point prediction, and an inviscid flow calculation. The com- 
plete analysis was applied to several separated afterbodies. The analysis accurately predicts 
the pressures and drag as a function of Reynolds number, Mach number, and afterbody shape 
for the comparisons made with data from afterbody-sting models tested by NASA. Thus the 
analysis provides an engineering approach for predicting nozzle afterbody drags. The small 
discrepancies found between analysis and data most likely result from inaccuracies in the 
separation point prediction, not the recirculating flow model. 

The limited comparisons of predicted and measured pressures over afterbodies with a flowing 
jet are encouraging. The predictions were made with the inviscid plume shape fixed. These 
calculations should, however, be repeated with a scheme to iterate between the external flow 
and the plume. The change in plume shape caused by the high pressures near the afterbody- 
jet junction should improve agreement between analysis and data. 

In closing, the viscous flow subroutine present herein provides a useful engineering tool for 
predicting pressures and drags over separated afterbodies. The coupling of the viscous flow 
subroutine to an inviscid flow calculation is highly sensitive, and extreme bare should be 
used when exercising the analysis. 


Langley Research Center 

National Aeronautic and Space Administration 
Langley Station, Hampton, Va., 28 June 1977 
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Figure 1 Schematic of the Separated Afterbody Flow field 



Figure 2 Comparison of Measured and Predicted Separation Angles 
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Figure 4 Differential Control Volume Analytical Model 
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Figure 9 Comparison of Calculated Discriminating Streamline Shapes For Circular Arc Nozzle Config- 
uration With f}= 17° and d^jD = 0,51 at 0. 94 
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Figure 1 0 Comparison of Calculated Discriminating Streamline Shapes For Circular Arc Nozzle Config- 
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Figure 12 Comparison of Calculated Discriminating Streamline Shapes For Circular Arc Nozzle Config- 

uration With B= 13.8° and djJD = 0.51 at 0.94 
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Figure 13 Comparison of Calculated Discriminating Streamline Shapes For Circular Arc Nozzle Config 

uration With j3 = 11° and d^jD = 0 \ 61 atM^ 0.8 
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Figure 14 Comparison of Calculated Discriminating Streamline Shapes For GrcularArc Nozzle Config - 
uration With (3= 11° and dJD = 0.61 at 0.90 
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Figure 16 Skin Friction and Pressure Gradient Effects on the Calculated Discriminating Streamline For 
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Figure 1 7 Discriminating Streamline Shape Effects on Afterbody Pressures For Circular Arc Nozzle 

Configuration With 0 = 17° and d^/D = 0.51 at M 00= q 6 
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Figure 18 Comparison of Predicted Separation Angles 



Figure 1 9 Major Components of the Calculation Scheme 
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Figure 20 Discrete Control Volume For Jet Entrainment Region 
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Figure 22 Schematic of the Flow Recirculation Regions 
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Figure 26 Convergence of the Viscous Inviscid Flow Iteration Scheme For A Circular Arc Conical 
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Figure 30 Comparison of Measured and Predicted Afterbody Pressures For A Circular Arc Nozzle Con- 
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Figure 35 Effect of Separation Location on Predicted Pressures and Drag of Circular Arc Nozzle Con- 

figuration With P= 13.4° and djjD = 0.51 at M 00 = 0. 798 
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Figure 36 Effect of Jet Entrainment on Discriminating Streamline Shape of a Circular Arc A fterbody 

at M 00 = 0. 6 and NPR =2. 91 
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APPENDIX 


AN AFTERBODY VISCOUS FLOW PROGRAM 
USER S MANUAL 

1.0 SUMMARY 

VISCUS is a computer subroutine package for analyzing the boundary layer flow separation 
occurring on aircraft afterbodies. The package is programmed to interface with existing com- 
puter programs that calculate external inviscid flow. Because VISCUS is a subroutine package, 
it cannot be executed without a controlling program. Therefore, to facilitate checkout, a 
main program — called the Stand-Alone Program — is being provided along with the VISCUS 
package. 

This appendix describes the VISCUS program interface and all routines in the package and 
also the Stand-Alone Program. The subroutine package and the main program can be ob- 
tained from COSMIC, 1 12 Barrow Hall, University of Georgia, Athens, Georgia 30601. 

2.0 GENERAL DESCRIPTION 

There are four major calculations in the VISCUS package: 

1 . calculation of the boundary layer 

2. determination of the separation pressure and the separation point 

3. calculation of the discriminating streamline 

4. combining of all components into an aerodynamic contour. 

The VISCUS subroutine package uses the modified Reshoto-Tucker method with the Ludweig- 
Tillman skin-friction law to calculate the boundary layer. In the first iteration, the boundary 
layer is calculated over the original contour. Subsequent calculations use the last calculated 
separated contour. 

Four different separation pressures can be calculated in the program, based on either the 
control volume criterion, Goldschmied’s criterion, the modified Page’s criterion, or Strat- 
ford’s criterion. The user selects the criterion, and the program calculates the separation 
point. 

After the second iteration, the separation point locations are averaged. Additional restraints 
on the movement of the separation point have been included to improve convergence. 

The separation point is used as the starting point for calculating the discriminating stream- 
line. The program uses an axisymmetric control volume calculation that accounts for chang- 
ing pressure, skin friction, and shear stress. There is also a “jet entrainment” option and jet- 
entrainment effects are included from the nozzle exit if the option is specified. 
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The primary output of the VISCUS subroutine package is the aerodynamic contour. And it 
is calculated by adding the boundary layer displacement thickness to the original or, if se- 
paration occurs, to the separated contour. The aerodynamic contour is returned to the in- 
viscid program which uses the revised contour to recalculate the external inviscid solution. 
See flow diagrams pages 75 to 77. 

3.0 USE OF VISCOUS SUBROUTINE PACKAGE 
3.1 GENERAL DESCRIPTION 


The Viscous Subroutine Package consists of all routines necessary to analyze boundary 
layer flow separation in the region of aircraft afterbodies. The Package has been designed to 
interface with an existing computer program system that calculates external inviscid flow. 
The “interface” is a fortran call that is inserted into the inviscid program after the program 
has computed surface pressure coefficients (C p ). The Viscous Subroutine Package uses the 
input contour and computed C p to determine an aerodynamic contour that is used to calcu- 
late inviscid pressures for the next iteration. The procedure is outlined below. 



> 


2nd TO Nth 
ITERATION 
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The number of iterations required to converge to a satisfactory solution is a function of 
the inviscid program and the configuration being analyzed. Experience has shown that ten 
iterations are usually sufficient. 

The following section describes the interface in detail. 

3.2 INTERFACE - FORTRAN CALL 

The user must include a call to the subroutine VISCUS, as described below. VISCUS 
is the controlling routine in the Viscous Subroutine Package, calling other subroutines as 
needed. 

CALL VISCUS (INT, FLOT, XA, YIN, CP, RJ,UJ, FLOTO, RADO, AM, THR, DEL, H, CFA, 
DELI, HI, RET, TAW, PTPT, FNN) 


Input Arguments 


INT 

integer array used to input integer scalars. Each scalar is assigned an 
array location. 

INTO) 

number of points used to described contour (max. = 500). 

INT (2) 

array location to start boundary layer calculation. 

INT (3) 

array location to start search for separation. 

INT (4) 

array location to end search for separation 

INT (5) 

jet entrainment option, 

= 0; omit jet entrainment effects, 

= 1 ; include jet entrainment effects. 

INT (6) 

array location of the nozzle exit. 

INT (7) 

smoothing option, 

= 0 ; do not smooth aerodynamic contour, 
= 1 ; smooth aerodynamic contour. 

INT (8) 

separation pressure option 
= 0: control volume 
= 1 ; Goldschmied 


= 2; modified Page 
= 3 ; Stratford 

FLOT 

real array used to input real scalars. Each scalar is assigned an array 
location. 

FLOT (1) - 

scale option, 

= 12.; input in inches, 
= 1.; input in feet. 

FLOT (2) - 

thermodynamic wall option, 
= 0. ; adiabatic wall, 

= 1 . ; isothermal wall. 
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FLOT (3) 

" 

free stream total pressure. 

FLOT (4) 

- 

free stream total temperature. 

FLOT (5) 

- 

free stream Mach number. 

FLOT (6) 

- 

interation counter. 

FLOT (7) 

- 

ratio of specific heats. 

FLOT (8) 

- 

initial boundary layer incompressible shape factor. 

FLOT (9) 


initial momentum thickness. 

FLOT (10) 

- 

not used. 

FLOT (11) 

- 

not used. 

FLOT (12) 

- 

shear layer spreading rate constant (input if INT (5) = 1). 

FLOT (13) 

- 

not used. 

FLOT (14) 

- 

not used. 

FLOT (15) 

“ 

the axial location of the separation point. If blank or zero, program 
will calculate separation point. 

XA 


array containing the axial locations of contour. 

YIN 

- 

array containing the radial locations of contour. 

CP 

- 

computed pressure coefficients array. 

RJ 

- 

inviscid plume radius array (input if INT (5) = 1 ). 

UJ 

- 

jet velocity array (starting at nozzle exit) (input if INT (5) = 1). 

Output Arguments 


FLOTO 

- 

array containing scalar output items. 

FLOTO (1) 

- 

pressure drag. 

FLOTO (2) 

- 

separation pressure (control volume). 

FLOTO (3) 

- 

separation pressure (Goldschmied). 

FLOTO (4) 

- 

separation pressure (modified Page). 

FLOTO (5) 

- 

separation pressure (Stratford) 

FLOTO (6) 


separation angle, degrees. 

FLOTO (7) 

- 

axial location of separation point. 

RADO 

- 

aerodynamic contour array. 

AM 

- 

Mach number array. 

THR 

- 

momentum thickness of the boundary layer array. 

DEL 

- 

displacement thickness of boundary layer array. 

H 

- 

compressible shape factor of the boundary layer array. 
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CFA 

DELI 

HI 

RET 

TAW 

PTPT 

FNN 


skin friction coefficient array. 

thickness of the boundary layer array. 

incompressible shape factor of the boundary layer array. 

Reynolds number based on momentum thickness array, 
adiabatic wall temperature array. 

ratio of boundary-layer-total-pressure to free-stream-total-pressure array, 
exponential shape factor of boundary layer profile array. 


3.3 INTERFACE SAMPLE 

The printout below is for a simplistic, inviscid program intended only as an example of 
how the interface between the inviscid program and the Viscous Subroutine Package might 
be handled. 


INVISCID PROGRAM 


DIMENSION 

1 

2 

3 

4 


INTI 8)* 
• RJ(500>, 

, THR ( 500 ) , 
, HI <500), 
* VAERO ( 500 ) 


PLOT (15)* 
UJ( 500) • 
DEL( 500), 
RETI500), 


X A ( 500 ) , 
FLOTO! 7), 
H ( 5 00 ) , 
TAW (500 ) , 


YINI500), 
RAD01500), 
CF M 500), 
PTPT (500), 


CP ( 500 ) 
AM( 500) 
DELI ( 500) 
FNN ( 500 ) 


READ (5, 
t 500 FORMAT ( 
WRITE! 
t»00 FORMAT ( 

1 

2 

3 

4 

READ (5 
510 FORMAT ( 
WRITE! 6 
fclO FORMAT! 

1 

2 

3 

4 

5 

6 

7 

8 


5001 HN,NMIN«NMAX* IJET *NEXT 
515) 

400) NN,NM1N,NMAX, I JET, NEXT 
NUMBER OF INPUT POINTS 
LOCATION TO START SEPARATION SEARCH 
LOCATION TO END SEPARATION SEARCH 
JET EN1RAINMENT OPTION 
LOCATION OF NOZZLE EXIT 
510) Z ,TWW,PT,TT,AMIN,GA,HIX,THRR,XSEP 
8F10.0) 

610) Z,TWW,PT,TT,AMIN,GA,HlX,THRR,XSEP 
SCALE OPTION 

THERMODYNAMIC WALL OPTION 
FREE STREAM TOTAL PRESSURE 
FREE STREAM TOTAL TEMPERATURE 
FREE STREAM MACH NUMBER 
RATIO OF SPECIFIC HEATS 
INITIAL BOUNDARY LAYER SHAPE FACTOR 
INITIAL MOMENTUM THICKNESS 
INPUT SEPARATION POINT (OPTIONAL) 


= • t 
=•, 
= • . 


15, / 
15, / 
15, / 
15, / 
15 ) 

,C,AlTER 

,Q,A1TER 

= t , 

= • 

= • 
s» 

=• 

=» 

= • 

=• 

= • 


F12.5, 
F12.5, 
F12.5, 
F12.5, 
F12.5, 
FI 2 « 5, 
F12.5, 
F12.5, 
F12.6, 


/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 
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9 • SHEAR LAYER SPREADING RATE = •, F12.5, / 

A • NO. OF ITERATIONS = •, F12.5, ) 

READ (5.520) I XAl I ) ,YIN( I ) ,1=1 ,NN ) 

520 FORMAT (2F10.0) 

WRITE (6* 620) ( XA( I ) * YINI!) ,1=1*NN ) 
o20 FORMAT I • X INPUT Y INPUT •, / (2FI2.5)) 


ANA IS THE ITERATION COUNTER 

ANA = 1. 

FIRST ITERATION 

DO 10 1=1, NN 
10 YAERO( I ) = Ylfil I ) 


20 CONTINUE 

CALL ROUTINE 10 CALCULATE 1NVISC1D SOLUTION 
CALL INV1S ( X A , YAERO , . . . , CP) 


SET UP ARRAYS FOR VISCOUS PACKAGE 


INK 1) 


NN 

I NT ( 2 > 


1 

I NT ( 3 ) 


NM1N 

INK A) 

= 

NMAX 

I NT ( 5 ) 

St 

I JET 

1NT ( 6 ) 

- 

NEXT 

INK 7) 


1 

INK 8) 


0 

F LOT ( 1 ) 


Z 

plot i 2 ) 


TWw 

PLOT (3) 

St 

PT 

PLOT ( A ) 

= 

TT 

FL0TI5) 


AMIN 

PLOT (6) 


ANA 

FLOT (7) 

= 

GA 

PLOT (6) 

zz 

HIX 

FLOT ( 9) 

— 

THRR 

FLOT ( 12) 

= 

c 

FLOT( 15) 

St 

XSEP 
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REMEMBER TO INPUT ORIGINAL CONTOUR TO V1SCUS 

CALL VISCUS I INT, PLOT, XA,YIN v CP»kJ»UJ*FLOTO'RAOO,AMtTHR, DEL 
1 ,H, CPA, DELI ,H1,RET,TAW,PTPT,FNN) 

CHECK IF ITERATION COMPLETE 

IF (ANA »GE» AlTER ) GO TO 999 

USE AERODYNAMIC RADII FOR NEXT INVISCID CALCULATION 

DO 30 1=1 ,NN 
30 YAEROI I ) = RADU(I) 

ANA = ANA+1. 

RETURN TO INVISCID SOLUTION 

GO TO 20 
* 

999 STOP 
END 


4.0 USE OF STAND-ALONE PROGRAM 


4.1 DESCRIPTION 

The stand-alone program consists of a main routine and the viscous subroutine package des- 
cribed in Section 3. The main routine reads input, calls the subroutine package, and prints 
output. The functions of the stand alone program are to define the input to the viscous pack- 
age, to provide an example of an interface with the viscous package, and to print output from 
one iteration of the viscous package. 


4.2 INPUT 
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DATA INPUT TO DECK B851 


DATE: 


CARD 


. TITLE: Input to Stand Alone Program 


FORMAT ( Namelist input 


) 


GENERAL INSTRUCTIONS : Follow general rules forlnameli st input. |The name of the input is NAMK, 


i 2 3 4 3 6 ? ■ » 10 It iz *3 •« 13 l« (7 It 19 20 21 2 2 23 24 2 3 2*27 2* Z* 30 3l 32 33 34 33 34 37 3* 39 40 41 *2 43 44 43 4« 47 4» 43 30 31 92 53 54 39 3* 37 36 39 SO *1 62 $3 «4 >3 tf 47 68 69 TO 71 72 73 7473 TS T T T* *3 


COLUMNS 

VARIABLE 

TYPE 

FIELD NAME 

REMARKS AND/OR DESCRIPTIONS 










Integer 

NN 

Number of input contour points fmax. = SO(Yl 






■ 

NAZ 

Array location to start boundary calculation. 






Integer 


Array location of minimum CP. The search for separation starts at 




this location. 






BSHB 

NMAX 

Array location of end of afterbody. The search for separation ends 


■HE 


at this location. 






Integer 

NEXT 

Array location of nozzle exit. 







IJET 

Jet entrainment option 




= 0 : no iet entrainment 




=/= 0 : iet entrainment 






Integer 

ISMOO 

Smoothing option. 

— 



= 0; no smoothing 




= 1 : smooth aerodynamic contour _ _ _ __ 






Integer 

1PRESS 

Separation pressure option 




~ = 0; control volume 




= 1 ; Goldschmied 




= 2; modified Page 




= 3: Stratford 






Real 

z 

Scaling option to convert input to feet units. 


— r* 


= 12. 9 innut inches 




= 1. . input feet 
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DATE 


DATA INPUT TO DECK 

CARD > TITLE: 

FORMAT ( ) 


GENERAL 

INSTRUCTIONS 










1 2 1 4 s • r a 

10 II IE 13 14 IS IS 17 

IS 11 10 El 22 23 2* 2S2S27ZB21 SO 31 32 33 3« 3S 3« 37 31 31 40 41 42 43 4443 4S 47 4S 41 SO 51 S2 S3 54 33 3S 37 3t 31 SO SI 12 S3S4S5 SSS7 SS S3 70 7i 72 73 74 73 7* t f r* *0 

COLUMNS 

VARIABLE 

TYPE 

FIELD NAME 

REMARKS AND/OR DESCRIPTIONS 


Real 

TWW 

Thermodynamic wall option 




= 0. , adiabatic wall 




= 1. , isothermal wall 






...Real.-. 

PT 

Free stream total pressure 





Real 

TT 

Free stream total temperature °R 





...Real ..... 

AMIN 

Free stream Mach number 





Real 

GA 

Ratio of specific heats 






Real 

HIX 

Initial boundary layer shape factor 






Real . 

THRR 

Initial boundary thickness. 


Real 

c 

Input only if jet entrainment option selected. Shear layer spreading rate 




constant. 


Real 

XSEP 

The axial location of the separation point. If blank or zero program will 




calculate the separation point. 






























1 




i 
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DATA INPUT TO DECK rrm 


DATE: 


CARD 


• TITLE: Input to Stand Alone Program 


FORMAT ( 
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DATA INPUT' TO DECK rxsI DATE: 

___________ 5 TITLE: Input to Stand Alone Program 


FORMAT ( ) 



7 /IE I ABLE 


COLUMNS TYPE FIELD NAME REMARKS AND/OR DESCRIPTIONS 
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5.0 DESCRIPTION OF ROUTINE 

This section briefly decribes all routines in the Viscous Subroutine Package. The routines 
are described in the order (top to bottom, left to right) outlined in the block diagram shown 
on the next page. 

VISCUS 

The VISCUS routine is the interface with the inviscid program and the controlling routine 
of the viscous package. The routine, in addition to these functions, calculates the total ve- 
locity from the inviscid surface pressure coefficient, calculates the cross-sectional area of the 
separated contour, and determines the axial location of the separation point. 

SHAPE! 

The SHAPEJ subroutine interfaces between VISCUS and SUMA and also sets boundary con- 
ditions for SUMA. 

SUMA 

The first derivative of the cross-sectional area versus axial location is calculated in SUMA 
routine. The routine fits a parabolic arc using a least-square fit before the derivatives are 
calculated. 

NEWBL 

The NEWBL routine calls the routine that calculates the boundary layer. After the boundary 
layer is calculated, the routine computes the component of the displacement normal to the 
x axis. 


BLC 

The BLC routine uses the modified Reshotko-Tucker method and the Ludweig-Tillman skin- 
friction law to calculate the boundary layer growth. The original contour is used for the first 
iteration. The separated contour is input for the other iterations. 

FIX 

The location of the minimum pressure coefficient is calculated in the FIX routine. The 
search region should be limited to the afterbody. 

POWER 

The axial coordinates are raised to integer powers and saved in a common statement. They 
are used in the separation pressure calculation. 
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SEPA 


The SEPA routine calculates the separation pressure. Four different methods are used. The 
“user” can select Stratford’s, Goldschmied’s, modified Page’s, or the Control Volume Criterion. 

PRFL 

The SEPA routine calls the PRFL routine to determine boundary properties based on the 
law-of-the-wall - the law-of-the-wake. 

FLUX 

The FLUX routine calculates integral boundary layer properties such as displacement and 
momentum thickness. 

INTEG 

The FLUX routine calls the INTEG routine to do the necessary numerical integration. The 
routine uses Simpson’s rule. 

SEP 


The SEP routine controls the calculation of the discriminating streamline and the aerodyna- 
mic contour. Before separation and after reattachment, the aerodynamic contour is the origi- 
nal contour plus the boundary layer displacement thickness. In the separated region, the 
aerodynamic contour is determined by adding the displacement thickness to the discrimin- 
ating streamline. 

B834 

The discriminating streamline is calculated in the B834 routine. A control volume technique 
is used. 

JET 

The JET routine calculates the effect of jet entrainment on the discriminating streamline. 
The JET routine is called only if the jet option is turned on. 

NEWRAP 

Both the B834 and the JET routines use the NEWRAP routine. The subprogram uses the 
Newton-Raphson method to predict the zero point of functions. 

SMINT-SMOOTH-MSMTH-RSMTH 

The SMINT-SMOOTH-MSMTH-RSMTH routines set up the contour for the smoothing rou- 
tine. They set boundary conditions and divide the contour into equal increments. 

SMTH 

The SMTH routine smooths the contour, using terms from the Taylor series expansion for a 
polynomial. 
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Viscous Subroutine Block Diagram 







6.0 GENERAL FLOW DIAGRAM 




CALCULATE BOUNDARY LAYER 


DETERMINE THE MINIMUM CP 


NORMAL TO CONTOUR AND TAKE 


ON THE AFTERBODY. THE MINIMUM 


COMPONENT NORMAL TO X-AXIS 


CP IS USED IN THE SEPARATION 




PRESSURE CALCULATION. 


CALCULATE SEPARATION 
PROPERTIES INCLUDING SEPARATION 
PRESSURE. FOUR SEPARATION 
PRESSURES ARE CALCULATED, 
STRATFORDS, GOLDSCHM lED'S MOD- 
IFIED PAGE, AND CONTROL VOLUME- 
THE USER CAN SELECT THE 
DESIRED CRITERIA. 


CALCULATE THE SEPARATION 
POINT. 



GENERAL FLOW DIAGRAM 
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